Activated events in glasses: the structure of entropic droplets 
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Using an effective potential approach, we present a replica instanton theory for the dynamics of 
entropic droplets in glassy systems. Replica symmetry breaking in the droplet interface leads to a 
length scale dependent reduction of the droplet surface tension and changes the character of the 
dynamical heterogeneity and activated dynamics in glasses. 
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At it's most basic level the puzzle presented by the 
existence of glasses is "how can any disordered arrange- 
ment of atoms be stable?" . Disorder implies myriads of 
structural possibilities - why don't the unguided thermal 
motions of the atoms in a glass take advantage of these 
possibilities and lead to ceaseless change and flow? Ac- 
tually, glasses do flow, but only very slowly. Thus our 
naive question must take on a more quantitative form- 
"why do materials with apparently high entropy density 
flow as slowly as they do?" Molecules in ordinary glass 
undergo displacements greater than an angstrom only on 
time scales greater than hours. Answering this quanti- 
tative question is tantamount to a theory of the glass 
transition. 

The most successful treatment of this issue to date is 
provided by the random first order transition (RFOT) 
theory of supercooled liquids and glassesQ, 0, IE 0, IE IE 
7| ■ Mean field calculations on these systems suggest a dis- 
continuous transition without latent heat at a Kauzmann 
temperature Tk, as found in infinite range spin glasses 
that lack reflection symmetry |8| , which is called "random 
first order" . The mean field approaches include the use of 
replicas[llllEIEIllEI3 and dynamical theories HE E3 
and exhibit 1-step replica symmetry breaking at Tk and 
a dynamical crossover temperature, Ta > Tk, for the 
onset of activated dynamics. Although a system with 
infinite range forces can break ergodicity while retaining 
a macroscopic configurational entropy, this is impossible 
for a system with finite range interactions 0,^3- Ergod- 
icity will be restored by "entropic droplets" , in which a 
region of the frozen solid is replaced by any of an ex- 
ponentially large number of alternatives. The entropy 
due to these many alternatives acts as a driving force for 
structure change which is, however, off-set by the free en- 
ergy cost of matching two alternative structures at their 
boundaries. This conflict gives the free energy barrier for 
slow activated motions. 

Many quantitative aspects of the glassy dynamics of 
molecular liquids have been predicted using the RFOT 
theory and are found to be in good agreement with ex- 
periment. The predictions include the relation between 
both the mean activation barriers in the liquid and glass 



regimes [l^ and the non-exponentiality of relaxation^) 
with the configurational heat capacity, the size scale of 
the mosaic structures ^j), which sho w up as "dynami- 
cal heterogeneities" in the laboratory [15| , as well as a 
quantitative prediction of the anomalous low tempera- 
ture properties of glasses, including their linear specific 
heat [la]- These successes of the RFOT theory encourage 
us to put the approximations presently used on firmer for- 
mal ground. This will allow the extension of the RFOT 
theory to systems more exotic than structural glasses 
such as colloidal systems having stripes or lamellae, i.e. 
micro-emulsions, which are predicted to be in the same 
universality class |l7|. or to glassy phases with strongly 
quantum degrees of freedom such as the electronic mo- 
tions of correlated electron systems |l8|. 

In this paper, using replica based methods, we study 
droplet dynamics in a frozen glassy system and show how 
it differs fundamentally from nucleation in ordinary first 
order transitions. These calculations yield a temperature 
and length scale dependent reduction of the droplet sur- 
face tension that strongly affects the dynamical hetero- 
geneity and activated dynamics of glassy systems. Specif- 
ically, we show that the transition states in a statistical 
field theory are instantons possesing additional replica 
symmetry breaking (RSB) in the entropic droplet and ar- 
gue that this new RSB state arises from configurational 
entropy fluctuations. Instanton like theories that ignore 
this complexity, particularly at the droplet interface have 
been presented earlier |3) UJl 120) l21| • Our results demon- 
strate that those approaches contain instabilities while 
the present RSB solution strongly supports the wetting 
picture put forward in[j|. 

We consider a system of interacting particles that be- 
comes glassy at low temperatures. Recent progress in 
the replica description of glass-forming systems with en- 
tropy crisis 0, IE E| allows us to formulate the RFOT 
theory in terms of a spatially varying collective variable, 
q a b (r) , that corresponds to the local Debye Waller fac- 
tor of molecular motions for liquid systems. q a t = 
in the ergodic liquid while q a b — > 0(1) for a ^ b in 
case of perfect long time correlations. Within replica 
mean field theory IE IE| the long time density corre- 



lations approach qpoS(k), with the equilibrium liquid 
structure factor being S (k), the mean particle density 
Po and the off-diagonal element q of q ab - Below the dy- 
namic transition temperature Ta, q becomes nonzero for 
the first time and an extensive configurational entropy 
S c occurs, that vanishes linearly at the Kauzmann tem- 
perature Tk < Ta- A heterogeneity of the glass, i.e. the 
emergence of correlated long time dynamics over certain 
length scales, is reflected in an r-dependence of q a b(r)- 

The mean field value and fluctuations of q ab are de- 
termined by an effective Hamiltonian, H[q\. In order 
to keep our calculation transparent we will not specifi- 
cally determine H [q] from an explicit calculation for su- 
percooled liquids (for an example see Ref.0), but start 
from a simple polynomial model in the same universality 
class 0: 

H = J feab] labqbcQca J (1) 

with h[q ab ] = \ (Vq ab f + \q\ h - ^q\ b + \q%. In 
Eq^ length scales are measured in units of the interpar- 
ticle distance and energy scales in units of a scale that 
determines the absolute value of the configurational en- 
tropy The parameters t, w > and y > 0, which are 
in principle all temperature dependent, become dimen- 
sionless. We take the primary T-dependence to reside in 
the quadratic term, t = T T r T ° where T determines the 
temperature scale where slow dynamics of the model sets 
in. We use both the absolute and reduced temperature 
variables, T and t, i.e. Tk or = t (Tk) etc. 

Our aim is to perform an instanton calculation to de- 
scribe the escape from a given metastable frozen solid 
configuration. This requires a technique that agrees 
in the homogeneous limit with other mean field cal- 
culations within the RFOT universality class fli| . 
but allows us to study the behavior for arbitrary 
values of the overlap between two states. This is 
achieved by the quenched effective potential approach 
of Franz and Parisi 6J . One considers a fixed overlap 
Pc(r) — pi(r)p2{r) between two particle density config- 
urations and determines the constrained partition func- 
tion Z(p c ,p 2 ) = (5(p c (r) - Pi(r)p2(r))) P - L by thermally 
averaging over one configuration. The effective potential, 
fi[p c ] = ~ T"(l°g Z(Pc P2))p 2 , is the average of the free en- 
ergy over the second configuration [(J and can be written 
as: 

[p c ] = 2 J d d rh \p c ] + F [ Pc ] . (2) 

Here F [p c ] = -T ^ / Dgexp (-/3H Pc [q])\ m ^ is the 
quenched free energy average for a system with Hamilto- 
nian, Eq^ in an external field p 2 c (r): H Pc [q] = H [q] — 
Sab / d d r P c (r) q a b (r). Replica theory gives the statisti- 
cal interpretation of p c and q ab . p c refers to the similarity 
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FIG. 1: Effective potential as a function of the overlap, p c for 
different temperature; inset: region of p c values where RSB 
occurs together with q*(t). 

with the initial state while q a b measures the configura- 
tional similarity between different states, regardless of 
the initial state. 

Before we perform the instanton calculation we sum- 
marize the homogeneous mean field theory of EqQ] fn 
Fig.l we show 51 [p c ] for r-independent p c where q ab in 
F [p c ] is determined via saddle point method. The min- 
ima of O [p c ] determine the mean field value of the Debye 
Waller factor q* = w +^^- A ± anc ) the configurational en- 
tropy fi [q*] = TS C @. We find t A = £ and t K = for 
the reduced dynamic crossover and Kauzmann temper- 
atures, respectively. The configurational entropy den- 
sity, s c = S c /V, vanishes linearly at the Kauzmann 
temperature, s c (t ~ tK) — j^Q* 2 (tic) (t — t K ), while 
s c (*a) = (tA)tA- For p c = q* , q ab is replica 

symmetric. The lowest replicon eigenvalue, A, is positive 

1/2 

for tK < t < tA and vanishes at t A like A oc (t A —t) . 
Since replica symmetry in the effective potential formal- 
ism corresponds to one step RSB in the conventional 
replica approach |22j, all these results are in agreement 
with the mean field theory of a RFOT. For p c values be- 
tween and a* an additional replica symmetry breaking 
in q ab occurs |J[23j (see inset of Fig.l). 

The il (p c ) of Fig.l clarifies the fact that between Tk 
and Ta the frozen solid with p c = q* is metastable with 
respect to the liquid (p c — 0) and higher in free energy by 
TS C . The system will try to reduce its energy by realizing 
this entropy via droplet nucleation. We argue that the 
transition state for this nucleation is determined by an 
instanton solution for the overlap profile p c (r), i.e. via 
j£H = with barrier Q (p c ) -Ct(q*) ; see also Ref-HU 

Assuming for the moment replica symmetry for the 
instanton, q ab (r) = (I — 5 ab ) q (r) we find the solution 
p c (r) = q(r), determined by Ap c (r) = ^jfj with 
V ( P c) = \v\ ~ f Pi + While V(p c ) and Q (p c ) dif- 
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for in detail, they have rather similar shape and agree 
for p c — q* . If replica symmetry holds, the problem 
then becomes almost identical to ordinary homogeneous 
nucleation|24j with transition temperature Tk and spin- 
odal Ta- Close to Tk a thin interface and a diverging 
droplet size result. The competition between surface and 
bulk terms determines the barrier 



AF = 4ir<T R 2 



An, 



-Ts c R J 



(3) 



The sole distinction from ordinary nucleation is that the 
driving force for the decay is entropic. For the surface 

1/2 



tension, the model EqQ] yields Co 
terface width is of order Iq 

radius is R ~ Iq t ^f - . The time scale for activated dy 



q* 2 {t K )t%, thein- 
2 J / 2 t K 1/2 and the droplet 
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at T, 
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namics therefore diverges like log — oc 
for spatial dimension d > 1. This result was first found 
in Ref.0 and also given in [l^ . l2li. Close to Ta no sharp 
droplet interface can be defined [2 5|. A dimensional anal- 
ysis of the instanton equation near Ta yields the length 
£4 ~ {pA {tA — t)) 1 ^ 4 , identical to the length scale in 
mode coupling theory above the upper critical dimension 
due = 6 26], see also Ref.|2f. Recently, using the effective 
potential approach, all these results for replica symmetric 
q a b were independently obtained by Franz |2l|. 

The self-generated randomness of a glassy system how- 
ever should modify the nucleation dynamics. A crucial 
question is whether the replica symmetry for q a t, is in- 
deed stable for the droplet. A stability analysis shows 
that the replicon eigenvalues, A, of the droplet solution 
are determined by 



(— V 2 + U (r)) tp (r) = Xt/j(r) 



(4) 



With- 



with r-dependent potential U (r) = d 2 

q 3=9* 

out the gradient term we would recover the result for A 
discussed above and replica symmetry is stable below Ta- 
The stability changes if we include, however, the gradi- 
ent term. Plainly, the replicon eigenvalues are identical 
to the eigenvalues of the usual fluctuation determinant in 
the decay of a metastable state. Thus, there is one nega- 
tive eigenvalue for Tk < T < Ta- Close to Tk it follows 
that A ~ — R~ 2 with droplet radius R. The unstable 
eigenvector is in the replicon direction and this instabil- 
ity occurs in addition to the expected instability with 
respect to the growth mode, i.e. the motion of the in- 
terface of the droplet. We further infer that the replicon 
eigenfunction, tp (r), is localized on the droplet interface. 

The instability of the replica symmetric solution entails 
an additional replica symmetry breaking as the overlap 
q c goes "over the hill", i.e. close to the interface. Allow- 
ing more general solutions for q a t, (r) — > q(r, x), where 
x G [0, 1] is the usual RSB parameter, we obtain from 

^y] = an d the saddle point equation for q a b, two 



coupled instanton equations: 

dv (p c ; 



V 2 Pc (r) 



and 



V 2 q(r, x) 



dp c (x) 



dV{q) 



Pc ( r ) + Pc (r) / q (r, x) dx (5) 



(1 - x) q 2 (r,x) -p 2 c (r) 



+ 



dq (r, x) 

dx q 2 (r, x) + 2q (r, x) I q(r,x'). 



(6) 



We have numerically solved these equations under the 
assumption of spherical symmetry and with the approx- 
imation that V 2 ~ li, justified for large droplet radius, 
i.e. for T close to Tk- We find that new solutions with 
one step replica symmetry breaking exist! In Fig. 2 we 
show our results for p c (r) and for q (r, x) below (qo (r)) 
and above ( gi (r)) the breakpoint, x. We find x to be 
r- independent. Outside the droplet, p c and the induced 
density correlations q (r, x) all approach the mean field 
value q* . However, inside of it, q{r,x) has a nontriv- 
ial distribution function. With weight x, q(r, x) is small 
(qo < p c ) and the transition state structures are configu- 
rationally distinct. However, with smaller weight, 1 — x, 
a few energetically well matching states contribute to 
making the transition. This demonstrates that dynam- 
ical heterogeneity in glasses differs fundamentally from 
ordinary nucleation. This instability originates in the 
interface, where the gradients are nonzero, but affects, 
via boundary conditions, the overlap distribution in the 
entire droplet. As T — > Tk we find x — > 1. At the 
Kauzmann temperature, where the replicon eigenvalue 
vanishes, the system is marginal in a new, heterogeneous 
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FIG. 2: Instanton profile for t = tn + 0.08(U - t K ) and 
y — w = 1 with breakpoint x = 0.975 and go,i = q(x ^x,r). 
Inset: surface tension as function of the droplet radius for the 
RS and RSB instantons. 



If we insert our solution into f2 [p c ] we obtain the in- 
teresting result that, compared to the replica symmetric 
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case, the bulk driving force for nucleation is unchanged. 
The new solution does however change the value of the 
surface tension. Close to Tk our numerical results are 
very well described by a (t) = CTo ^1 — b ^-t K ) w ^ n l ar S e 
values for the constant b (b ~ 10 for w = y = 1 and b ~ 5 
for u> = 3, y = 1). The reduction of the surface tension, 
shown in the inset of Fig. 2 as function of droplet radius 
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Ts m , is caused by the fact that RSB recognizes the 
existence of structures whose interface matches are less 
costly. 

The driving force for the decay of the metastable frozen 
solid is entropic. As usual for a replica based theory, 
S c is only the mean value of the logarithm of the num- 
ber of states. While fluctuations of the entropy density 
vanish for infinite systems, (5s c ) 2 — c c /V , they become 
relevant in finite subsystems like small droplets. For ex- 
ample we find for the configurational specific heat ca- 
pacity at Tk'- c c — 2q 3 T^ . Configurational entropy 
fluctuations, i.e. fluctuations of the driving force for 
nucleation, act similarly as do field fluctuations in the 
random field Ising model (RFIM). In the RFIM simple 
statistical arguments yield a reduction of the surface ten- 
sion a (R) = (To (l — const. i?~ 3 / 2 ). The effect is caused 
by fluctuations, corresponding to a wetting of the inter- 
face by states that make the matching between inside 
and outside the drop easier, a (R) computed this way 
behaves very similar to our numerical results presented 
above. In the context of the RFIM, Villainjij} proposed 
that one can consider such a a (R) as the first step of 
a renormalization group argument to the surface tension 
and obtained a cx R ~ from the integration of the 
flow equation. For an RFOT this leads to a droplet 
radius R* cx (i — tx )^ 2 ^ d and a Vogel-Fulcher law: 
t = tq exp [DTk/(T — Tic)] for the activated dynamics 
as advocated in Ref . Our results strongly support this 
wetting picture in entropic droplets. For practical appli- 
cations in liquids, where it is hard to approach Tk very 
closely, the present first step of the full renormalization 
group calculation, may well be sufficient. 

In summary we have shown that the transition states 
responsible for the slow activated dynamics in glassy sys- 
tems are given by instantons with new replica symmetry 
breaking (RSB). This new RSB is due to a changed dis- 
tribution of states at an interface between a frozen solid 
and a state which explores various metastable states. We 
note, the present treatment also allows naturally the in- 
clusion of effects that will lead to barrier height fluctua- 
tions and non-exponential glassy relaxation. 
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